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Ch ' Determining ground states of correlated electron systems is fundamental to understanding novel 

O' 

, ^, ' phenomena in condensed matter physics. A difficulty, however, arises in a geometrically frustrated 

system in which the incompatibility between the global topology of an underlying lattice and 

> 

lO ' local spin interactions gives rise to macroscopically degenerate ground states , potentially prompt- 

(N' 

^O ■ ing the emergence of quantum spin states, such as resonating valence bond (RVB) and valence 

cn 

f~^ . bond solid (VBS). Although theoretically proposed to exist in a kagome lattice ~ one of the most 

o, 

C^ , highly frustrated lattices in two dimensions (2D) being comprised of corner-sharing triangles - such 

quantum-fluctuation-induced states have not been observed experimentally. Here we report the 



^ ■ first realization of the "pinwheel" VBS ground state in the S = ^ deformed kagome lattice antifer- 

romagnet Rb2Cu3SnFi2. In this system, a lattice distortion breaks the translational symmetry of 
the ideal kagome lattice and stabilizes the VBS state. 
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Theoretically, significant progress has been made toward understanding the ground state 
of the quantum [S = |) kagome lattice antiferromagnet via various approaches. Diverse 
classes of novel states have been proposed; those include a gapless U(l)-Dirac-spin-liquid 
state^'^, a gapped-spin-liquid state^"^, and a VBS state^"— . One leading approach is built 
upon a quantum dimer model, which was first introduced to describe spin dynamics of 
singlet pairs (dimers) in high-T^ superconductors^^'^-^. Based on this model, two kinds of 
quantum states promptly emerge, namely the VBS and RVB states^i^i^i""— . The VBS state 
has long-range dimer-dimer correlations that break the translational symmetry of the ideal 
kagome lattice with spin-1 singlet-to-triplet excitations, whereas the gapped RVB state, in 
which different dimer coverings resonate among themselves, has short-range dimer-dimer 
correlations and topological order with continuum spin-i excitations (spinous). Since these 
quantum states are non-magnetic and lack static order, one needs to look for their distinct 
features in magnetic excitations by means of inelastic neutron scattering to distinguish them 
in a real system. 

Experimentally, realizations of these states have been extensively studied in quasi-lD^^"— 
and Shastry-Sutherland 2D systems^^. In contrast, the rarity of model systems has hitherto 
precluded experimental attempts to observe these states in the kagome lattice. In recent 
years, the discoveries of the kagome lattice in herbertsmithitei^"— and distorted kagome 
lattice in volborthite^ii^I have been generating a lot of excitement and debate over their 
plausible ground states. However, the lack of single crystals and Zn-Cu intersite disorder 
in herbertsmithite hinder the study of these systems. In this work, we present a single- 
crystal study of a stoichiometrically pure S = ^ deformed kagome lattice antiferromagnet 
Rb2Cu3SnFi2 with full occupancy of spin-i at the Cu sites. Combined neutron scattering 
measurements on a large single crystal and advanced numerical analysis provide the first 
direct evidence for the pinwheel VBS state. 

The magnetic S = ^ Cu^"*" ions of Rb2Cu3SnFi2 reside at the corners of the triangles of 
the deformed single- layer kagome lattice, and are inside distorted octahedral cages formed 
by F^ ions. At a room temperature, the system crystallizes in the hexagonal R3 space- 
group with the lattice constants a = 13.917(2) A and c = 20.356(3) A^^ (Figure la). The 
slightly distorted kagome lattice yields four antiferromagnetic exchange interactions (Fig- 
ure lb), which are characterized by different Cu^^— F^— Cu^+ bond angles ranging from 
124° for the weakest bond to 138° for the strongest bond. The system does not exhibit 



magnetic order down to 1.3 K— i^, which is confirmed by our elastic neutron scattering on 
a powder sample. At low temperatures, magnetic susceptibility rapidly decreases toward 
zero, a signature of a non-magnetic, spin singlet {S = 0) ground state. Furthermore, low 
temperature magnetization sharply rises when the applied field H exceeds a critical field 
He ~ 20 T, indicative of singlet-pair breaking and closing of a singlet-to-triplet gap. Finite 
magnetic susceptibility observed at the base temperature when the field is applied perpen- 
dicular to the crystallographic c-axis is attributed to the mixing of the singlet and triplet 
states via the antisymmetric Dzyaloshinskii-Moriya (DM) interaction^^. Therefore, to a first 
approximation, the spin Hamiltonian is given by: 

n = J2 [JiA ■ ^3 + ^ii ■ Si X s,] , (1) 

nn 

where Jij are the nearest-neighbor interactions for four different bonds Ji > J2 > J3 > J4 as 
labelled in Figure lb, and Dij = Di, D2, -D3, and D^ are the DM vectors, whose in-plane 
components point toward the center of a triangle and out-of-plane components alternately 
point into and out of the kagome plane^^. The spatial distribution of the non-uniform 
exchange interactions suggests the stabilization of the pinwheel VBS state^- in which singlet 
dimers are formed between two spins bonded by Ji ~ 20 meV^^. This pinwheel VBS state 
is different from the VBS state theoretically proposed for the ideal kagome lattice, since 
it is stabilized by the lattice distortion and the long-range dimer-dimer correlations do not 
break the translational symmetry of the lattice. A theoretical study of the J4-depleted 
kagome lattice^, which is defined by Ji = J2 = J3 = J and J^ = aJ, shows that the most 
energetically-favored VBS state proposed for the ideal kagome lattice, which corresponds 
to a = 1, is destabilized against this pinwheel state even at a ~ 0.97, highlighting the 
importance of geometrical frustration which leads to the diversity and closeness of exotic 
phases in the quantum kagome lattice. 

To verify the pinwheel VBS state, we performed inelastic neutron scattering on a single 
crystal sample (see Methods). The scattering intensity can be described by the following 
equation: 

/(Q,M= (^)V(Qm^,q;Qm), (2) 

where (^) = 72.65 x 10^^ barn/yU^, /(Q) is a magnetic form factor for Cu^"*", S(hu, q; Q^) 
is a dynamical structure factor, and Q = Qm+Q- An energy scan at the Brillouin zone center 
Qm = (0, 2, 0) (Figure 2d) clearly shows two singlet-to-triplet gaps. The convoluted fit with 



the resolution function yields the gap energies Ai = 2.35(7) meV and A2 = 7.3(3) meV. 
The thermal renormalization of the spectral weight and damping of the spin excitations 
(Figure 3) are typical of triplet modes in an interacting dimer systent^I. In addition, energy 
scans in an applied magnetic field (Figure 4) reveal the Zeeman splitting of the lower gap 
and no splitting of the upper gap. These results demonstrate that, whereas the upper gap is 
a result of the excitations from the singlet to S^^^ = triplet states, the lower gap, which is 
two- fold degenerate at zero field, is due to the excitations from the singlet to 5*^^^^ = ±1 triplet 
states (S'fot represents a triplet spin state of a single dimer). The Zeeman splitting of the 
S^g^ = ±1 triplet gap can be described by the following relation Ai^±{H) = Ai — g±fiBS^otH, 
where /i^ is the Bohr magneton. For S^^^ = 1 and —1, we obtained g^ = 2.58(30) and 
g_ = 2.34(20), respectively. Extrapolating the S^^^ = 1 line to intercept the H-a.xis, i.e. 
Ai +(// = He) = 0, yields the critical field He ~ 21 T, consistent with the magnetization 
datai2^. The magnitude of the DM vector can be estimated from Moriya's calculation^^ as 
\D\/Jav ~ Ag/g ~ 0.2 (Jav = \ J2i ^i)- Both Lande gf-factor and \D\ are in good agreement 
with the values obtained from fitting the susceptibility data to the exact diagonalization for 
the 12-site kagome cluster— 1^. 

For Rb2Cu3SnFi2, six dimers in the unit cell give rise to six triplet modes. In the limit 
J2 — Ja ^ 0, these modes can be characterized by irreducible representations of the cyclic 
group C3 as four E and two A localized (flat) modes (see Figure 6 in ref. 26 and Figure S3 in 
the Supplemental Information). The observed triplet excitations belong to the lowest-energy 
E mode. The partial lifting of this E-mode triplet degeneracy can be attributed to the DM 
interaction. For a single dimer, one can easily show that the DM interaction raises the energy 
of the S'^.t = state, while retaining the degeneracy of the two S'^^^ = ±1 states. To explain 
the observed large bandwidth and small S^^^^ = ±1 triplet gap, we performed the dimer series 
expansion on the pinwheel VBS ground state using the model Hamiltonian given in equation 
(1). We assume that the strength of the DM vectors scales with the respective exchange 
interactions, i.e. dz = Dt/Ji {dp = D^/Ji) where i =1, 2, 3, and 4 for the out-of-plane 
(in-plane) component. In our first-order dimer expansion (see Figure S3), the out-of-plane 
component dz of the DM interaction lowers the S^^^ = ±1 branch at F point, which strongly 
suggests that the observed two branches are the lowest-energy E-mode in the S^^^ = ±1 and 
sectors. On the other hand, the in-plane component dp only determines detailed structure 
of further energy splitting, and thus is omitted in the higher-order calculations. The linked 



cluster expansion algorithn>2^ was used to obtain a longer series. The lowest-energy triplet 
excitation spectra were then calculated using the Dlog-Pade approximation up to eighth 
order in the inter-dimer exchange and DM interactions. Under the constraints Ji > J2 > 
J3 > J4 and a fixed Jav, we obtain the best fit to the excitation spectra (Figures 2a-c) 
for Ji = 18.6 meV, J2 = 0.95Ji, J3 = 0.85Ji, J4 = 0.55Ji, and d^ = 0.18. The agreement 
between the measured and calculated triplet dispersions verifies the pinwheel VBS state. The 
integrated dynamic structure factor (Figure 2e) was calculated to fifth order, assuming equal 
contributions from two domains, one domain with pinwheel arms in clockwise direction and 
the other with pinwheel arms in the counterclockwise direction. These two pinwheel motifs 
correspond to two equivalent structural domains, which are a mirror image of each other 
(see Figure S5 in Supplemental Information). The intensity pattern with the periodicity four 
times the reciprocal lattice vector correlates with the intra-dimer distance, which is one forth 
of the in-plane lattice constant. The maximum scattering intensity around (2, 0, 0), (0, 2, 0) 
and (2, 2, 0) also supports the pinwheel VBS state in Rb2Cu3SnFi2 (see Supplemental 
Information). 

In contrast to triplet excitations in other dimer systems, the highly-dispersive triplet 
bands and the small triplet gap (Ai/Ji ~ 0.13) in Rb2Cu3SnFi2 are due to the strong inter- 
dimer interactions, the DM interaction, and the unique pinwheel arrangement of dimers. 
The bandwidth of the S^^^ = branch is denominated by J2 — J^-, therefore, the rather small 
Ji ~ ^2/2 is essential to explain its relatively large bandwidth. It is also worth noting that 
the S'fot = branch is virtually unaffected by dz and has a gap minimum at K point, as 
predicted by the bond operator mean field theory without the DM interaction^^. On the 
other hand, dz changes the gap minimum of the S^^^ = ±1 branch from K point to F point, 
and significantly lowers the gap energy. If the strength of dz exceeded 0.25, the gap at F 
point would be closed and magnetic order would instead be observed. 

The macroscopically degenerate ground states in the kagome lattice are particularly sen- 
sitive to small perturbations, and thereby could inherently give rise to a number of different 
states in a real system. In Rb2Cu3SnFi2, our results reveal the pinwheel VBS state that 
has been suggested to be in close proximity to other phasesiii^^. Recently, we have studied 
other related compounds, in which a different ground state was observed^^. In particular, 
Cs2Cu3SnFi2 exhibits magnetic order but possesses the low-temperature phase that shares 
the crystal structure with Rb2Cu3SnFi2, which suggests similar magnetic interactions. The 



magnetic order in Cs2Cu3SnFi2 could, therefore, be ascribed to (i) a small difference in 
exchange coupling and/or (ii) a large DM interaction^. In this case, the former could play 
a more prominent role in determining the ground state since the DM interactions in these 
systems are essentially comparable^^. In addition, other phases could be obtained by tuning 
external parameters such as magnetic field and pressure. Further experimental and theoret- 
ical studies are desirable to determine precisely the driving mechanism that leads to these 
different ground states. A more comprehensive understanding of this mechanism could shed 
some light on many remaining questions regarding the quantum kagome lattice, particularly 
in relation to the sought-after quantum spin liquid. 

Methods 

Experiments 

Inelastic neutron scattering experiments were preformed on the HER and GPTAS triple-axis 
spectrometers operated by the Institute for Solid State Physics, University of Tokyo. Vertically 
focusing pyrolytic graphite (PG) crystals were used to monochromate the incident neutron beam. 
For the measurements on the HER spectrometer, the scattered neutrons with a fixed final energy 
of 5 meV were analyzed by the central three blades of a seven-blade doubly focused PG analyzer. 
A cooled Be or oriented-PG-crystal filter was placed in the incident beam, and a room-temperature 
Be filter in the scattered beam to remove higher-order contamination. For the measurements on 
the GPTAS spectrometer, vertically focused (horizontally flat) PG crystals were used to analyze 
scattered neutrons with fixed final energies of 14.7 meV and 13.7 meV. The horizontal collimation 
sequence of 40' - 80'-sample-80' - 80' (LR) or 40' - 40'-sample-40' - 80' (HR) was employed 
with a PG filter placed in the scattered beam. A single crystal of mass 0.94 g was mounted in the 
hkO and hOl zones. The single crystal sample used in this study was synthesized using the method 
described in ref. 23. The sample was sealed in an aluminum container in a He gas environment 
for heat exchange, and cooled by a closed cycle ^He cryostat. For the measurements in a magnetic 
field, the sample was aligned in the hkO zone, and the field was applied along the crystallographic 
c-axis, perpendicular to the kagome plane, using a vertical magnet with a maximum field of 5 T. 
The sample was then cooled by a liquid ^He cryostat equipped with the magnet. 



Data Analysis 

In Figure 2a, the measurements were done around (2, 0, 0) for energies between 2 meV and 4.5 
meV, and around (2, 2, 0) between 4.5 meV and 9 meV. In Figure 2b, the measurements were done 
around (2, 2, 0) along two equivalent high symmetry directions; for energies between 2 meV and 4 
meV, the scans were taken along (H, H, 0), and between 4 meV and 9 meV along (H, 6-2H, 0). In 
Figure 2a, the triplet branch around (1.5, 0, 0) could be a result of higher-order corrections to the ex- 
change interactions due to small atomic displacements as evidenced by weak superlattice reflections; 
we have ruled out the contribution of other triplet modes (see Supplemental Information). All scans 



were analyzed using an empirical dispersion relation, hujaici) = \/ ^a ~ (~1)"(^q!,x^x +^a,«^«)) 
where a = 1 or 2 for two non-degenerate branches, and q = {kx,ky) is a wave vector away from 
a Brillouin zone center Qm- Within the resolution of the instrument of ~ 1 meV (GPTAS-LR), 
the out-of-plane dispersion is flat, attesting to the two-dimensionality of the system, and hence we 
set Voi,z = 0. The observed peaks were fit with narrow, resolution-limited Lorenztians, convoluted 
with the resolution function taking into account the empirical dispersion relation. The line shapes 
are governed by the convolution with the resolution function. The peak positions in a unit of A^^ 
were obtained by fitting Va and an offset Aq (see Figure S2 in Supplemental Information). 
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FIG. 1: Crystal structure of Rb2Cu3SnFi2. a A crystal structure in the ab plane shows the 
connectivity of the Cu^"*" ions (red) forming a deformed kagome lattice. In b, dimers (blue bonds) 
form the pinwheel VBS state. The exchange interactions are labelled as Ji > J2 > J3 > J4. 
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FIG. 2: Energy-momentum contour maps. Intensity contour maps were measured at T = 3 K 
around (2, 0, 0) and (2, 2, 0) along two high symmetry directions, a along k^ and b along ky. c 
Triplet dispersions were measured along the path shown in the inset. Solid lines represent the best 
fits to the dimer series expansion described in the text. Error bars reflect the differences of various 
Fade approximations, d An energy scan at Q^ = (0, 2, 0) (F— point) and T = 5 K shows two spin 
gaps at 2.35(7) meV and 7.3(3) meV. e The calculated dynamic structure factor integrated over 
triplet excitations shows high scattering intensity around (2, 0, 0), (0, 2, 0) and (2, 2, 0), consistent 
with the experimental data. 
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FIG. 3: Representative energy scans and temperature dependence, a-c Energy scans of the lower 
gap at three temperatures. Sohd lines denote the fits to the dispersion relation described in Meth- 
ods, convoluted with the instrumental resolution function. Peak widths are resolution-limited at 
the base temperature and become broader with increasing temperatures. Temperature dependence 
of the two triplet gaps is shown in d. 
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FIG. 4: Triplet excitations in a magnetic field, a-c Energy scans measured at T = 5 K show the 
splitting of the lower gap at three different fields. Solid lines denote the fits to the dispersion relation 
described in Methods, convoluted with the instrumental resolution function. Field dependence of 
the triplet gaps is shown in d. 
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Supplemental Information 

Superlattice reflections 

Below Tg ~ 215 K, X-ray diffraction stiows superlattice reffections at positions displaced from 
the fundamental Bragg peaks by commensurate wave vectors (±1/2, ±1/2, 0), (±1/2,^1,0), and 
(±1,=F1/2,0) (Figure Sl(a)). The largest intensity of the superlattice peaks is roughly 50 times 
smaller than the fundamental Bragg peaks, indicative of small atomic displacements. The absence 
of an anomaly in magnetic susceptibility at Tg, as well as the relatively large superlattice intensity 
compared to that in other inorganic spin Peierls systems such as CuGeOs and TiOCl, suggests 
that this structural phase transition is not the spin Peierls transition. Hysteretic behavior of 
the superlattice intensity measured by neutron scattering (Figure Sl(b)) suggests a first order 
structural transition at Tg. Below Tg, the hexagonal lattice is retained but the R3 space-group 
changes to P3 or P3 with an enlarged unit cell. 

Triplet excitations around (1.5, 0, 0) and (2.5, 0, 0) 

Dispersive branches of the triplet modes were also observed around (1.5, 0, 0) and (2.5, 0, 0) 
(see also Figure 2a). A long constant-energy scan around (2, 0, 0) (Figure Sl(c)) shows triplet 
excitations around (1.5, 0, 0) and (2.5, 0, 0), which disappear at high temperatures. Figure 
Sl(d) shows the singlet-to-triplet gap at (1.5, 0, 0) with Ai = 2.21(9) meV, consistent with the 
triplet gap at (2, 0, 0) within the experimental uncertainties. The presence of these extra triplet 
branches and the spin gap at (1.5, 0, 0) indicates the change in the periodicity of the magnetic 
exchange interactions due to the enlarged unit cell of the low-temperature phase as evidenced by 
the superlattice reflections. The small atomic displacements at low temperatures slightly alter the 
exchange-interaction network, and thereby give rise to a new periodicity of the triplet modes. The 
scattering intensity of the triplet modes around (1.5, 0, 0) and (2.5, 0, 0) is about five times smaller 
than that around (2, 0, 0). The fit in Figure Sl(c) already takes in account the magnetic structure 
factor of the Cu^+ ion. In order to account for this half periodicity, one needs to perform the dimer 
series expansion on 24 dimers or 4 pinwheels in the enlarged unit cell, which is beyond the scope 
of this study. However, this correction is small, and does not change the pinwheel arrangement of 
dimers. The origin of this weak branch will be further discussed in the next section. 
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Calculations of triplet excitation energies at M point 

The possibility that the observed dispersive branches around the M points (1.5, 0, 0) and (2.5, 
0, 0) are due to other triplet modes is ruled out as described below. 

We calculated the whole triplet excitation energies at M point using the dinier series expansion 
method up to eighth order. For E modes, we obtained reasonable convergence of the Dlog-Pade 
approximations. The two lowest energies of the E modes with S^^^ = O(ibl) at M point are 5.5 and 
6.5 meV (5.1 and 7.3 meV), which are much higher than 2.4 meV. For A modes, a bare energy series 
indicates the tendency of softening at M point as a function of A, which denotes a scaling parameter 
for the inter-dimer and DM interactions, i.e. Jj = a^A Ji where i = 2, 3, 4 and Qj's are constants. In 
order to check if it is a second order phase transition, we examined the Dlog-Pade approximations 
of the eighth order series, in which we found that most of the Dlog-Pade approximations have 
poles and the critical exponent for the vanishing gap is estimated as ~ 0.2. Since this critical 
exponent is much smaller than the known gap exponent of 0.71 for two-dimensional systems, we 
conclude that the softening found in the eighth order series of the A mode energies is artificial 
due to the shortness of the series. We have no reliable energy estimations of the A modes because 
of the aforementioned artificial softening. However, our structure factor calculation indicates that 
scattering intensities of the A modes are much smaller than those of the E modes, which rules out 
the possibility that the scattering intensity at M point in the constant energy scan with hw ~ 2.4 
meV comes from the A modes. Therefore, we conclude that the extra mode at (1.5, 0, 0) and (2.5, 
0, 0) is due to the magnetic super-structure, which corresponds to the enlarged structural unit cell 
evidenced by the superlattice refiections. 

Representative constant-energy scans and dispersion relations (Figure S2) 

Constant-energy scans were performed around (2, 2, 0) at T = 3 K with (a) hw = 3.5 meV along 
ky, (b) hoj = 3.0 meV along kx, and (c) fuxj = 6.5 meV along kx- The solid lines show the fits to 
the dispersion relation described in Method, convoluted with the instrumental resolution function, 
(d) and (e) show dispersions of the singlet-to-triplet excitations along kx and ky, respectively. In 
(d), the HER data are shown in green diamonds. The dotted lines represent the model dispersion 
described in the text with a universal Va^i = 35 meV/Afor a = 1, 2 and i = x,y. 
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Characteristic features of the energy spectra 

We here present the calculated energy spectra within the first-order perturbation, for comparison 
with the bond-boson calculation employed in ref. 26. This lowest-order perturbation calculation 
provides comprehensive explanations for some characteristic features of the excitation spectra of 
the pinwheel VBS state. 

We first focus on the exchange interaction, omitting the DM interaction, and describe how the 
dispersive modes emerge. As pointed out in ref. 26, when J2 = J4,, a pinwheel, hatched regions in 
Figure lb, becomes an orthogonal dimer ring, which prohibits a triplet dimer from propagating. 
Hence, the spectra are non-dispersive (flat) as shown by black solid (E-mode) and dashed (A-mode) 
lines in Figure S3(a), since a triplet dimer is confined around a triangle formed by J3. The energies 
of E- and A- modes are given by Ji — f and Ji + ^, respectively. A non-zero value of J2 — J4 leads 
to the breakdown of this orthogonal nature and the triplet-excitation spectra become dispersive as 
shown by red curves in Figures S3 (a). 

We now turn to the effect of the DM interaction. We first consider only the out-of-plane 
component dz, where S^^f is a conserved quantity. Within the first-order perturbation, the effective 
Hamiltonians for the Sf^^ = ±1 sectors depend on d^ but that for the Sf^^ = sector does not. 
The calculated spectra for dz = 0.09 are shown in Figure S3(b). The S^^^ = 1 branches (blue 
curves) are affected by dz 7^ 0, whereas the Sf^^ = branches (red curves) are unchanged. Figure 
S3(c) shows the result of dz = 0.18, in which the two lowest-energy branches denoted by thick lines 
resemble the measured dispersion curves. We finally include the in-plane component dp] the results 
for dz = dp = 0.18 are shown by green lines in Figure S3(c). We found that dp only determines 
detailed structure of the spectra and has a very small effect on the two lowest-energy branches. 

Dynamic structure factor 

The dynamic structure factor S{q,uj) is defined by 



^.> (3) 



where N denotes the total number of spins, rj^^ the position of the ^^^ spin on the i^^ unit 
cell, and ipg the ground state. Introducing a complete set of momentum eigenstates |V'( y )) of 
the Hamiltonian, we can rewrite the dynamic structure factor as S{q,cj) = X^a'^aC^''^) with 
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S'a( g ,a;) = lA{q)S{uj — EA{q) + Eg). In order to obtain the perturbative estimation of SaCq ,(^), 
we first utilize the hnked cluster expansion algorithm to calculate series expansions of the exclusive 
matrix elements (see ref. 29) and the effective Hamiltonian Heff- The energy E^i^) is obtained 
as an eigenvalue of Heff, and the weight IaCq) is calculated from the exclusive matrix elements 
and an eigenvector of H^ff- 

In Figure 2(e), we show the fifth order result of the integrated structure factor of triplet dimer 
(TD) states, ItdCq) = "l^AeTD I ^^'^^C^ ^^) ~ l^AgTD ^a( y )• Since we carry out numerical 
diagonalization of the fifth-order effective Hamiltonian to get eigenvectors required in the weight 
calculation, our calculated IaCq) is not in the form of a power series and thus we do not per- 
form series analysis such as Fade approximation. Although our weight calculation contains no 
extrapolation scheme, it can reproduce the experimental feature of high scattering intensity in the 
Brillouin zones centered at (2, 0, 0), (0, 2, 0), (2, 2, 0) and so on, consistent with the experimental 
result. 

Now we describe how the spectral weight Itd is allocated to individual triplet-dimer modes. 
Figures S4(a) and (b) show the calculated structure factor S{q, oj) along two equivalent high sym- 
metry directions (kx), (H, 0, 0) and (H, 2, 0), respectively. Similarly, Figures S4(c) and (d) show 
the calculated structure factor along the other two equivalent high symmetry directions (ky), (H, 
H, 0) and (H, 6-2H, 0), respectively. These high symmetry directions correspond to the measuring 
paths shown in Figures 2a and b. We find four branches in the measuring energy range {fuv < 9 
meV); we call them, from the bottom, first, second, third, and fourth branches. The first and 
fourth (second and third) branches belong to Sz = ±1(0) mode, and are denoted by dashed (dot- 
ted) lines. In Figure 2a, the measurements were done for 2 meV < hco < 4.5 meV along the (H, 
0, 0) direction and for 4.5 meV < fkj < 9 meV along the (H, 2, 0) direction. In Figure 2b, the 
measurements were done for 2 meV < fkj < 4 meV along the (H, H, 0) direction and for 4 meV 
< fiuj < 9 meV along the (H, 6-2H, 0) direction. These measuring ranges are indicated by white 
rectangles in Figure S4. 

In Figures S4(a) and (b), the presence of the intense first branch is in agreement with the 
experiment. However, the weak second branch and intense third branch in (b) are inconsistent 
with the experiment. Note that in the dispersion fitting we have assigned the observed higher- 
energy mode to the second branch. Furthermore, while the second branch appears more intense 
than the third branch in (a), it becomes less intense in (b). (Similar inconsistencies were also 
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observed along ky in Figures S4(c) and (d), though they are less obvious). These results imply 
that (2, 0, 0) and (2, 2, 0) are singular points of S{q,u}). This comes from the fact that the 
second and third branches are degenerate at T point. For dz = 0, all four branches converge at F 
point. The degeneracy of the first and forth branches at F point is, however, lifted by the non- 
zero dz, whereas that of the second and third branches remains. Hence, we infer that our model 
Hamiltonian, in which we neglect the in-plane DM interaction and the larger-structure of exchange 
interactions, is not sufficient for the purpose of discussing the scattering intensity of an individual 
mode. Since the in-plane DM interaction does not resolve the degeneracy at F point (see Figure 
S3), the incorporation of the larger-structure of exchange interactions is likely required to settle 
the intensity issue of the second and third branches. On the other hand, the out-of-plane DM 
interaction dz = 0.18 determines the wave functions of the first and fourth branches, for which 
we have no singular behavior around F point and the larger-structure of exchange interactions is 
expected to be irrelevant. We note that in Figure S4(b) the fourth branch is as intense as the first 
branch, which also contradicts the experimental result. This calculated result of the fourth branch 
with rather strong intensity may suggest that the in-plane DM interaction, whose strength dp is 
expected to be the same order as dz, strongly affects the triplet-dimer wave functions. (Note that 
our omission of the in-plane DM interaction is based on the consideration of excitation energies.) 
In summary, there is still an inconsistency between the experimental data and numerical calcu- 
lations. We believe that this inconsistency is not essential, but is artifact due to the degeneracy at 
F point. Additional terms in the spin Hamiltonian, which resolve the degeneracy of the second and 
third branches at F point, could remove the singularity of S{q,oj) and alter the dynamic structure 
factors of the second and third branches. Further studies which include such terms, as well as 
the incorporation of an adequate extrapolation into our perturbation calculation, are apparently 
necessary to settle this issue. 

Magnetic domains: clockwise and anticlockwise pinwheel states 

We believe that in Rb2Cu3SnFi2 there are two magnetic domains which give rise two different 
arrangements of dimers forming clockwise and anticlockwise pinwheel states as shown in Figure S5. 
The presence of these two magnetic domains, which are a mirror image of each other, corresponds to 
two structural domains. These two different arrangements of dimers generate a distinct scattering 
pattern. Figure S6 shows the integrated dynamic structure factor of the clockwise pinwheel state, 
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which after combined with that of the anticlockwise state gives the scattering pattern shown in 
Figure 2e. 
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Figure SI (a) X-ray photographs show the fundamental Bragg peaks and the superlattice peaks 
indicated by the arrows at positions (±1/2, ±1/2,0) away from the fundamental peaks, (b) 
Intensity of the superlattice peak (5, 3.5, 0) measured with neutron scattering on a single crystal 
shows the hysteretic behavior at Tg. (c) A constant-energy scan &t hw = 2.4 meV reveals the 
triplet excitations around (1.5, 0, 0), (2, 0, 0) and (2.5, 0, 0). (d) An energy scan at (1.5, 0, 0) 
yields Ai = 2.21(9) meV. Lines are fits to the empirical dispersion relation described in Methods 
convoluted with the resolution function. 
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Figure S2 (a)-(c) Representative constant-energy scans along k^ and ky. Lines are fits to the 
empirical dispersion relation described in Methods convoluted with the resolution function. 
Dotted lines in (d) and (e) show the empirical dispersion relation with a universal Va = 35 
meV/A. 
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Figure S3 First-order results of the triplet-dimer excitation spectra. For the right axis, we use 
Ji = 18.6 meV. 
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Figure S4 Calculated dynamic structure factor S{q, uj) for the energy range of ^ < 9 meV, where 
the weight /a is estimated by the fifth-order calculation and energy E\ by the eighth-order series 
with the Dlog-Pade analysis. White rectangles denote the measuring ranges plotted in Figure 2a 
and b. 
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Figure S5 Two magnetic domains, (a) a clockwise pinwheel state and (b) an anticlockwise 
pinwheel state. 



Figure S6 The integrated dynamic structure factor is calculated for the clockwise pinwheel state. 
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